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In terms of spin coherent states we have investigated topological defects 

in 2D S = 1 (pseudo)spin quantum system with the bilinear and biquadratic 
isotropic exchange in the continuum limit. The proper Hamiltonian of the 
model can be written as bilinear in the generators of SU(3) group (Gell-Mann 
matrices). The knowledge of such group structure allows us to obtain some 
new exact analytical results. Analysing the proper classical model we arrive at 
different skyrmionic solutions with finite energy and the spatial distribution 
of spin-dipole and/or spin-quadrupole moments termed as dipole, quadrupole, 
and dipole-quadrupole skyrmions, respectively. Among the latter we would 
like note the in-plane vortices with the in-plane distribution of spin moment, 
varying spin length, and the non-trivial distribution of spin-quadrupole mo- 
ments. 

PACS numbers: 71.10.-w; 75.10.Hk; 75.10.Jm 



1 



I. INTRODUCTION 



Different topological defects play an important role both in low-energy (spin excitations, 
domain walls, superfluidity/superconductivity) and high-energy physics from heavy ion col- 
lisions to cosmological scenariosiHl. Theoretical approach to its description traditional for 
strongly correlated systems like quantum (pseudo)spin ones, starts from either (pseudo)spin 
Hamiltonian with subsequent reduction to either classical models with solutions like in-plane 
or out-of-plane vortices, and skyrmions. The latter represent the solutions of non-linear o- 
model with classical 2D Hamiltonian 
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H n = J d 2 r 



Hi) 2 



i=l 



(1) 



for the vector field n(r) = {sin 6 cos $, sin 6 sin cos obtained by Belavin and Polyakov@ 
more than two decades ago. A renewed interest to these unconventional spin textures is 
stimulated by high-T c problem in doped quasi-2D-cuprates and quantum Hall effect. 

The skyrmion spin texture consists of a vortex-like arrangement of the in-plane compo- 
nents of spin with the ^-component reversed in the centre of the skyrmion and gradually 
increasing to match the homogeneous background at infinity. The spin distribution within 
classical skyrmion is given as follows 

r 2q - \ 2q 

* = W cosg= r2g + A v ( 2 ) 

or for q = 1 



2rA 2rA r 2 - ,\ 



In terms of the stereographic variables the skyrmion with radius A and phase (po centered 
at a point z is identified with spin distribution w(z) = jz^, where z = x + iy = re %Lp is 
a point in the complex plane, A = Xe l6 , and characterized by three modes: translational, 
or positional Zo-mode, "rotational" ^-mode and " dilatational" A-mode. Each of them corre- 
sponds to a certain symmetry of the classical skyrmion configuration. For example, #-mode 
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corresponds to a combination of rotational symmetry and internal U(l) transformation. 
Classical skyrmionic energy E q = 8nqJS 2 is proportional to its topological charge and does 
not depend on its radius. 

Other well known solutions of isotropic and anisotropic 2D Heisenberg model are the 
in-plane and out-of-plane vortices@~i which have the energy logarithmically dependent on 
the size of the system. The in-plane vortex is described by the formulas $ = qip, cos 8 = 0. 
The 6(r) dependence for the out-of-plane vortex cannot be found analytically. 

Non-linear cr-model can be addressed as classical continuum limit of 2D Heisenberg fer- 
romagnet with isotropic spin-Hamiltonian 

H=-jJ2^S i+ b. (4) 

i,b 

The simplest quantum generalization of skyrmionic solutions could be obtained in frames 
of spin coherent statesi, where the wave function of the quantum spin system, which maxi- 
mally corresponds to classical skyrmion, is a product of spin coherent states. In the case of 
spin s = \ 

* sfc (0) = n[ cos | e ^ It) +sm|e-^ (5) 

i 

where 9i = arccos[(r 2 — A 2 )/(r 2 + A 2 )]. The coherent state implies a maximal equivalence 
to the classical one with the minimal uncertainty of spin components. Actually, every 
on-site spin in a lattice is assumed to be subjected to a molecular field H(r) oc n(r) = 
{sin#cos$, sin#sin$, cos#} which spatial distribution forms a skyrmionic texture. 

The coherent state approach appears to be rather simple for the s = 1/2 spin systems. 
Indeed, on the one hand, spin-Hamiltonian for s = 1/2 quantum system is restricted to have 
isotropic, or anisotropic bilinear Heisenberg exchange form like (^). On the other hand, the 
trial wave function (|5|) is simply parameterized by the vector field n(r). Some quasiparticle 
properties of quantized skyrmion in the s — 1/2 model are addressed irM 

The situation becomes more involved for the S > 1 (pseudo)spin systems, where, gen- 
erally speaking, we have to deal with additional non- Heisenberg terms in (pseudo) spin- 
Hamiltonian and several vector fields to parameterize the trial wave function like (||). A 
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principal difference between the S = ~ and S > 1 quantum systems lies in what con- 
cerns the order parameters. The only single-site order parameter in the former case is 
an average spin (dipole) moment (S XjVjZ ), whereas in the latter one has additional "spin- 
multipole" parameters like " spin-quadrupole" (spin nematic) averages ({SiSj}), where 
{SiSj} = SiSj + SjSi. Hence, we may expect in S = 1 quantum spin systems differ- 
ent topological defects with spatial distribution of not only nonzero spin (dipole) moment 
(dipole skyrmions), or spin-quadrupole moment (quadrupole skyrmions), but more involved 
dipole-quadrupole skyrmions with a nonzero distribution of both order parameters. 

Interestingly, that in a sense, the s = 1/2 quantum spin system is closer to the classical 
one (S — > oo) also characterized by the single-site vector order parameter than, for instance, 
the S = 1 quantum spin system with its eight site order parameters. In a whole, we should 
expect for the S > 1 (pseudo)spin systems an appearance both of unconventional topology 
and the complicated order parameter textures^. 

The main interest in S = 1 quantum spin systems is provoked by the S = 1 quantum spin 
chains displaying the Haldane gapS. Several special examples of the 2D S = 1 spin systems 
have been extensively discussed earlier in connection with the study of the static and dynamic 
properties of anisotropic Heisenberg 2D magnets with a single-ion anisotropy (see e.g.0E!). 
Some aspects of the topological structure of vortices in 2D S = 1 systems were discussed by 
different authors in connection with 3 He problemlli'0 and triplet superconductivity^. The 
quantum "spin nematic" phase of the fully isotropic 5=1 system was addressed recently 
by Ivanov and KolezhukS. 

Our interest in this field was motivated by a problem of a description of soliton-like 
excitations in quasi-2D cuprates in frames of a novel scenario proposed by one of authors0. 
The model considers doped cuprates as a system of singlet local bosons moving in a lattice 
formed by hole centers CuO\~ . Such a center has a complex ground state multiplet 1 Ai g — 1 E u 
which can be described by a pseudo-spin S — 1. 

In this paper we make use of the spin coherent (SC-) state approach to describe the 
skyrmion-like topological defects in S = 1 quantum (pseudo)spin 2D systems with isotropic 



non-Heisenberg spin-Hamiltonian. These solutions are a special case of so-called CP at spinors 
discussed in Refs.llllll, though its authors focused their interest only in CPi and CP3 models, 
rather than the CP2 spinors of our model. Our interest is mainly focused on the unconven- 
tional quadrupole and dipole-quadrupole skyrmion-like static solutions. 

In Sec. II we address the isotropic bilinear-biquadratic Hamiltonian for the S = 1 quan- 
tum (pseudo)spin systems, the parameterization of the trial wave function, the SU(3)-model 
approach, and the reduction procedure to the classical continuum limit of the S = 1 model. 
In Sec. Ill the unconventional skyrmion-like solutions are analyzed, including the known mag- 
netic (dipole) skyrmionS, and unusual quadrupole and dipole-quadrupole skyrmions with a 
non-trivial spatial distribution of dipole and/or quadrupole (pseudo)spin order parameters. 

II. CLASSICAL DESCRIPTION OF THE S = 1 QUANTUM (PSEUDO)SPIN 

SYSTEMS 

In general, isotropic non-Heisenberg spin-Hamiltonian for the 5 = 1 quantum 
(pseudo)spin systems should include both bilinear Heisenberg exchange term and biquadratic 
non-Heisenberg exchange term: 

H = —Ji SiSi +v — J2 ^^(SiSi +v ) 2 = (6) 

i,rj i,rj 
„ . 3 

Vre.ci fo.ol 

l+T)J 



—J\ SiSi +v — J2 ^^({SkSj}i{SkSj} 

i,r) i,r) k>j 



where Jj are the appropriate exchange integrals, J\ = J\ — J%/2, J 2 = t/ 2 /2, i and rj denote 
the summation over lattice sites and nearest neighbours, respectively. In our spin-1 model 
we use trial functions 

v>= n c ^>*= n ( a iti)+Mj))ipi m 

jdlattice jalattice 

Here j labels a lattice site and the spin functions ipi in cartesian basis are used: ip z = |10 > 
and ip X:V ~ (| 11 > ±|1 — 1 >)/\/2- The linear (dipole) spin-operator is represented by simple 



matrix: 

< ipi\Sj\i/} k >= -ie ijk , 
and for the order parameters one easily obtains: 

< S >= -2 [a, b], < {SiSj} >= 2(5ij - aiaj - bibj) (8) 

given the normalization constraint a 2 + b 2 = 1. Thus, for the case of spin-1 system the order 
parameters are determined by two classical vectors (two real components of one complex 
vector c — a + ib from ([/])). The two vectors are coupled, so the minimal number of dynamic 
variables describing the S = 1 spin system appears to be equal to four (see Refill and Sec. 2. 2 
below). Hereafter we would like to emphasize the director nature of the c vector field: ip{c) 
and ip(—c) describe the physically identical states. 

The dipole or magnetic skyrmions in the spin-1 systems were addressed in the paper0. 
The structure of the order parameter admits the existence of more general types of solutions 
which are purely quadrupole ("electric") or mixed dipole-quadrupole ("magneto-electric") 
ones. One should note that, in common, the length of the spin vector in S=l model must 
not be fixed. The order parameters structure is responsible for another important property 
of the S=l systems: it allows the existence of more than one topological quantum number. 

A. SU(3)-symmetry model: Gell-Mann operators and effective Hamiltonian for S=l 

model 

Three spin-linear (dipole) operators 61,2,3 an d five independent spin-quadrupole operators 
{Si, Sj} — 7^S Sij at S — 1 form eight Gell-Mann operators being the generators of the SU(3) 
group. Below we will make use of the appropriate Gell-Mann 3x3 matrices A^ k \ which 
differ from the conventional only by a renumeration: \^ = A^, A^ 2 -* = A^ 3 \ A® = A®, 

A (4) = A ( 5)) A (5) = _ A (2) > A (6) = A (4) ; A (7) = A (l) ; A (8) = A (7) _ Firgt three matrices A (l,2,3) 
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correspond to linear (dipole) spin operators: 

/0 \ / i\ 



A (1) = S x 



-i 
\0 i / 



A® = S„ 




\-i 0/ 



A (3) = S, 



/ -i \ 

i 
\0 0/ 



while other five matrices correspond to quadratic (quadrupole) spin operators: 
/0 0\ /0 1\ 



A (4) = 



1 

Vo i o/ 



;A® = -{S X S,} 





Vi o o/ 



/0 1 o\ 



A® = -{5^} 





/I 





o \ 


A (7) = -^ + ^-2^) = i= 





1 









Vo 





-2/ 



A (8) = ^ 2 - Si 



1 

Vo o o/ 

/l 0\ 

0-10; 
\0 0/ 



+ S y + S z 



IE 



with being a unit 3x3 matrix. 

The generalized spin-1 model can be described by the Hamiltonian bilinear on the SU(3)- 
srenerators 

^ = -EE^ (fc)A S- 

i,i+r) k=l 

Here i, rj denote lattice sites and nearest neighbors, respectively. This is a S = 1 counterpart 
of the S — 1/2 model Heisenberg Hamiltonian with three generators of the SU(2) group or 
Pauli matrices included instead of eight Gell-Mann matrices. 

In frames of isotropic bilinear-biquadratic model we study in this work the 8x8 matrix 
Jkm is assumed to be diagonal with elements Jn = J 2 2 = J33 = J\, and J 44 = J 55 = 
... = J$8 = Ji- Fully isotropic SU(3) model with J\ = J 2 corresponds to a ferromagnetic 
version of so-called Uimin-Lai-Sutherland mo It should be noted that the isotropic 

in a real space Hamiltonian with J\ ^ J2 can be considered as the anisotropic one in the 
8-dimensional SU(3) group space, and the symmetry of the model breaks to the subgroup 
SO(3)cSU(3). In other words, the breaking of the condition J\ = J%{— J) can be considered 
as an appearance of the exchange anisotropy in the 8-dimensional phase space. To describe 
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this anisotropy one might introduce the ratio A = Jij J\. Hereafter, we shall associate the 
limiting cases A = and A — > oo with purely magnetic and electric solutions, respectively. 
Such an effective anisotropy in S — 1 systems differs strongly from the real spatial exchange 
anisotropy in S — 1/2 systems which results in preferred spin orientations. 

If one considers the magnetic field parallel to z-axis or the anisotropy of exchange pa- 
rameters in a real space the symmetry breaks to SO (2) C SO (3). However, given the definite 
relations between the anisotropy constants and exchange integrals this model can be reduced 
to spin-1/2 isotropic model. This case merits the separate examination. 



B. The classical continuum limit of S=l model 

Having substituted our trial wave function (^) to (H) provided (S(1)S(2)) = 
(S(1)){S(2)) we arrive at the Hamiltonian of the isotropic classical spin-1 model in the 
continuum approximation as follows: 



H = J x d 2 f 



i=l 



J 2 d' 2 r 



4(J 2 - Ji) 



S\ 2 d 2 f, (9) 



where S = — 2[a*b] = (S). The effective exchange anisotropy defines the third "gradient less" 
term in the Hamiltonian that breaks the scaling invariance of the model. Such an effect in 
S — 1/2 system appears due to the real spatial exchange anisotropy which defines the 
magnetic length!. 

The spin-1 model differs from the spin-1/2 model due to the appearance of the additional 
("nonmagnetic") degrees of freedom. When (S) = (e.g., if one of a, b vectors turns into 
zero) we have a non-zero part of classical Hamiltonian (proportional to J 2 ) an d can get 
the nontrivial configurations of non-zero vector. We shall call this configuration "electric 
skyrmion". It should be described by one vector with the fixed length, so the topological 
classification of such solutions is completely analogous to that of the spin-1/2 classical solu- 
tions, which are described by the order parameter being a fixed-length spin vector. Below, 
we shall study this solution, deriving it via reducing our biquadratic model to the non-linear 



S 



0(3)-model. The topological charge of the classical electric skyrmion with 6 = can be 
defined by usual formula as follows 

Q = — / d 2 fe utl (a* [a v * a J) = — / rdrdtp(a* [a r * -aj) , (10) 

where the subscripts denote derivative. 

In the continuum limit for J\ = J 2 = J the Hamiltonian @ can be transformed into the 
classical Hamiltonian of the S'L r (3)-symmetric scale-invariant model: 
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H =\ J I ^E(v#< fc V) 2 ]. 
^ k=i 



where we make use of the single-site wave function in the form as follows: 

/i*iexp(i$i)\ 



if, 



R\ 2 = l, (11) 



i? 2 exp(i$ 2 ) 

with i? = {sin B cos r], sin B sin r], cos B}. In accordance with the director nature of the a, b 
vector fields we have to vary the angles Q,rj, $j in the range (0,7r). The Hamiltonian can 
be rewritten as follows: 

H lsotr = 2J J rf 2 f{(VB) 2 + sin 2 B(Vr/) 2 + 
+ sin 2 Bcos 2 B cos 2 r^(V^i) 2 + sin 2 r/(W 2 ) 2 + sin 4 6 cos 2 r] sin 2 ^(V^i - W 2 ) 2 } , (12) 

where we have introduced = $i — $3, \l/ 2 = $ 3 — $ 2 . For an isolated (pseudo)spin system 
the $3 phase is arbitrary, hence the minimal number of dynamic variables describing the 
S = 1 (pseudo)spin system equals to four (=450). However, for a more general situation, 
when the (pseudo)spin system represents only the part of the bigger system, and we are 
forced to consider the coupling with the additional degrees of freedom, the $ 3 phase turns 
into a non-trivial parameter: $ 3 = Q^ip. The topological solutions for our Hamiltonian (|l~2"|) 
can be classified at least by three topological quantum numbers (winding numbers): phases 
rj, \l/ li2 can change by 2ir after the passing around the center of the defect. It should be noted 
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that the director nature of the a, b vector fields implies the possibility for winding numbers 
to take half-integer values. The appropriate modes may have very complicated topological 
structure due to the possibility for one defect to have several different centers (while one 
of the phases 77, ^1,2,3 changes by 2tt given one turnover around one center {r\,(p\), other 
phases may pass around other centers (r^, (pi)). Each center in a multi-center defect can be 
considered as a quasiparticle. In this connection it should be noted that the spin-1 model 
differs from the spin-1/2 one by the fact that the former in common assumes quasiparticles 
of different types due to the existence of different topological quantum numbers for different 
centers. 

Finally, it should be noted that the above model approach can be extended to the S > 1 
isotropic (pseudo)spin systems. In the case of spin S one has to calculate averages like 
< Si 1 ..Si n > where n — 1, .., 2S. One can suggest that in our Hamiltonian the combinations 
Ci x ..Ci n with a fixed- length vector c should appear. Now let us notice that for any k e N the 
model with discrete Hamiltonian 



Hkk — —Jkk 



.9=1 



(13) 



in the continuum limit can be reduced to 



Hh — Jk 



d 2 f 



jWfc=i 



Jk- 



(14) 



Making use of V(/g) = fVg + gVf provided |c| =const, we come to 

3 



H k = kJ k \c\ 



2k~2 



St 



8=1 



(15) 



It is a non-linear 0(3)-model. Its skyrmionic solutions differ from the conventional ones by 
the energy, due to the term k\c\ 2k ~ 2 . In our spin-1 case k = 2, c = n. So, if the reduction 
of spin-S quantum model to the quasi-classical one ( |13"D will be made by means of some 
parameterization of the trial wavefunction it will be easy to obtain the skyrmionic solutions 
of the model with finite energy kJkq\c\ 2k . 
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III. UNCONVENTIONAL SKYRMIONS IN S = 1 (PSEUDO)SPIN SYSTEMS 



A. Dipole skyrmions 

One important case of the spin-1 model when J 2 = (purely Heisenberg Hamiltonian) 
also has skyrmionic solutions, which were found earlier in Ref.0. When the a, b vectors are 
perpendicular to each other (a _L b), the model also reduces to the nonlinear 0(3)-model. 
The solution from Ref.0 is described by the following formulas (in polar coordinates): 

\Pla = (e z sin 9 — e r cos 9) sin <p + e v cos if ; 

V2b = (e z sin 9 — e r cos 9) cos <p — e, ' fi sin (p . (16) 

The fixed-length spin vector is distributed in the same way as in the conventional skyrmion 
with topological charge q = 1 (@). However, unlike the usual skyrmions, the solutions (|16D 
have additional topological structure due to the existence of two vectors a and b. Going 
around the center of the defect the vectors a and b can make iV turns around the spin 
vector oc [a x b] . Thus, we can introduce two topological quantum numbers: N and g0 In 
addition, it should be noted that q number may be half-integer. Hereafter we shall call the 
skyrmionic solutions with the only non-zero magnetic component as the dipole or magnetic 
ones. 



B. Quadrupole skyrmions 

Magnetic skyrmions as the solutions of purely Heisenberg (pseudo)spin Hamiltonian0 
were obtained given the restriction a _L b and the lengths of these vectors were fixed. These 
restrictions lead to the pure magnetic solution and enabled to use a subgroup for the topo- 
logical classification^. It was SO(3) which is generated by A.^, matrices forming 
the Heisenberg bilinear term. 

Hereafter we address another situation with purely biquadratic (pseudo)spin Hamiltonian 
( Ji=0) and treat the non-magnetic (electric) degrees of freedom. The topological classifi- 
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cation of the purely electric solutions is simple because it is also based on the making use 

— * 

of subgroup instead of the full group. We address the solutions given a || b and the fixed 
lengths of the vectors, so we use for the classification the same subgroup as above. 
The biquadratic part of the Hamiltonian 

3 



Hk 



can be rewritten as follows 



d 2 f 



.i,j=l 



(17) 



Hi 



biq 



I 



Jo I d 2 f 



(18) 



where a = an, b = (3n, and a + i(3 = exp(m), k e R. We denote n = 
n{sin 9 cos sin 6 sin $, cos 9}. Using simple formula V(fg) = fVg + gV f together with 
the normalization constraint |n| 2 =const, we reduce the expression for H biq to the familiar 
nonlinear 0(3)-model: 



H, 



biq 



2J 2 \n\ 2 I d' z r 



]T(Vn, 



(19) 



Its solutions are skyrmions, but instead of the spin distribution in magnetic skyrmion we 
have solutions with zero spin, but the non-zero distribution of five spin-quadrupole moments 
^ A. (4,5,6,7,8)^ or ({SiSj}) which in turn are determined by the distribution of the a(n) vector: 

r 2 9 _ A 2 9 



$ = qip + <3> ; cos 9 



r 2q _|_ \2q 



(20) 



with a classical skyrmion energy 



Eei = 16nqJ 2 . 



(21) 



The distribution of the spin-quadrupole moments {{SiSj}) can be easily obtained: 

(r 2q + \ 2q ) 2 cos 2 + (r 2q - \ 2q ) 2 sin 2 qtp _ 



( r 2 9 + X 2qy 



{r- 2q + \ 2q ) 2 sin 2 qip + 



„2q 



X 2q ) 2 cos 2 qip 



{r 2q + A 2 ?) 2 
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2r 2 ^sin2gy 2 4r 2 ^ 



( r 2<j _|_ ^2 9 )2 ' \ 2/ ( r 2g + ^2 g )2' 

/ro C1X 2(A^-^)r^singy /ro C1X 2(A 2 « - r^A" cos q<p 

One should be emphasized that the distribution of five independent quadrupole order pa- 
rameters for the electric skyrmion are straightforwardly determined by a single vector field 
n(r). The phase factor a+i/3 can be arbitrary, because it is not included in the Hamiltonian. 
It may be written as follows: a + if} = exp(ipip). Such classes of solutions are physically 
equivalent in frames of our model. 

C. Dipole-quadrupole skyrmions 

In this subsection we would like to write out some solutions of fully SU(3)-symmetric 
isotropic model for (pseudo)spin Hamiltonian @ with J\ = J 2 . We do not derive all the 
possible solutions of this model because it merits to be a subject of the separate examination. 
Our main goal now is only to illustrate the richness of the model using as examples the 
simplest solutions. The solutions of this model may be classified taking into account whether 
they have each of the winding numbers to be zero or not. Below we will briefly consider 
two simplest classes of such solutions (the choice of classes is defined by the simplicity of 
integration) . 

First of all we shall consider simplest solutions of the equations minimizing energy func- 
tional. One type of skyrmions can be obtained given the trivial phases $12. If these are con- 



stant, the R vector distribution (see (ITU)) represents the skyrmion described by the usual for- 
mula @. All but one topological quantum numbers are zero for this class of solutions. It in- 
cludes both dipole and quadrupole solutions: depending on selected constant phases one can 
obtain both "electric" and different "magnetic" skyrmions. The substitution $x = <3> 2 = $3 
leads to the electric skyrmion which was obtained as a solution of more general SU(3)- 
anisotropic model in the previous section. Another example can be $1 = $2 = 0, $3 = tt/2. 
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This substitution implies b\\Oz, a\\Oxy, S\\Oxy, and S = sin cos ©{sin ,7, — cost], 0}. Nom- 
inally, this is the in-plane spin vortex with a varying length of the spin vector 

- _ 2r\\r 2 - A 2 | 
11 (r 2 + A 2 ) 2 ' 

which turns into zero at the circle r = A, at the center r = and at the infinity r — > oo, and 
has maxima at r = \{y/2 ± 1). In addition to the non-zero in-plane components of spin- 
dipole moment (S x>y ) (or (A^ 1 ' 2 ^)) this vortex is characterized by a non-zero distribution of 
(pseudo)spin-quadrupole moments (A^ 6 ' 7,8 )). 

Here we would like to emphasize the difference between spin-1/2 systems in which there 
are such the solutions as in-plane vortices with the energy having a well-known logarithmic 
dependence on the size of the system and fixed spin length, and spin-1 systems in which 
the in-plane vortices also can exist but they may have a finite energy and a varying spin 
length. The distribution of quadrupole components associated with in-plane spin-1 vortex 
is non-trivial. Such solutions can be named as "in-plane dipole-quadrupole skyrmions". 

Other type of solutions we get given the phases \I>i = Qi^p, ^ = Q2P with two integer 
winding numbers Qi,2 and r) = i](r),Q = 0(r). Easiest way to obtain the solutions of 
variational equations is to put one of these functions to be constant. This way we arrive at 
six types of solutions. The first one 6 = corresponds to the trivial spinor 



(23) 



Other solutions being written in complex form with z = re lip /\ (z = z*) are as follows: 

( 1 \ 

(it): 8 = it/2 + nk,k e Z,i;(z) = arctan[|z| ( " +,3) ]; V 





/i?iexp(i$i)\ 




(i) : V = 


R 2 exp(i$ 2 ) 


= 




\ Rs J 


W 



exp(iqi(p) 



(1 + | 2 |2(gi+g 2 ))i/2 



( — l\k^(qi+q 2 ) 
/ 

(24) 



{Hi) : 77 = irk, k G Z, Q{z) = arctanfl^l 91 ]; ip 



(l + |^|2?l)l/2 



/{-l) k z qi \ 


1 



(25) 
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(iv) : r\ = 7r/2 + nk, k G Z, 6(2) = arctanf^l 92 ]; ^ 



1 



(1 + | Z |2?l)l/5 



/ \ 

(-l)fe^a 

V 1 / 



; (26) 



(v) : 77 = vr/4 + tt/c/2, ie2,gi = -g 2 = 9, 6(2) = arctan[|z| 9 ]; ^ = - /2 



/^/V2\ 

V 1 / 

(27) 



Id 29 - 1 

(vi) : 77 = 7r/4 + irk/2, k G Z,q 1 = q 2 = q, Q(z) = arccos[ , [0 „ , ^ ]; ip 



z \2q + X J' i- n + UI2,) 



Vi - M 29 / 

(28) 



The energy for solution (vi) is .Eg — 167rgJ, while for others ((ii) — (v)) E}. = Aa^J 
with a 2 = qi + q 2 , a 3 = <?i, «4 = <?2 and 05 = g. Similar solutions, as so-called CPjv spinors 

/ a \ 



rj(z) 



1 



+ jVr 2 «)V 2 



(29) 



\ w 

(q is a winding number, a is any constant) are discussed in Ref. though its authors 
focused their interest only in CPi and CP3 models, rather than the CP2 spinors of our 
model. 

Above we address the dipole-quadrupole solutions with the only angular parameter hav- 
ing the z(r, <£>)-dependence. Hereafter, we consider the nontrivial solution with two winding 
numbers and unconventional z-dependence of angular parameters. To this end we make use 
of a -^-spinor in another form as follows: 
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/ sin(|) exp(if) (008(f) sin(f ) + i sin(f ) cos(f )) \ 
ip = sin(|) exp(i|)(sin(|) sin(f) - icos(f) cos(f)) 
V cos(f) / 

In the system where S z is diagonal it takes form 

/sm(|)cos(f)exp(^)\ 



(30) 



V 



/ 



sm(f)sin(f)exp(i2=2) 
cos(f) 

Such a parametrization is opportune when we are interested in solutions with mean-values 
of S z and Q zz independent of ip. The density of the energy functional with J\ — J 2 — 1 can 
be written out as follows: 



l r/ - 



sin 2 9 r 



{(Va) 2 + 2 sin $(Va, V/3) + (V/3) 2 }+ 



+ sin 2 -{(V$) 2 + cos 2 $ sin 2 -(V/?) 2 }]. 
2 2 



(31) 



The substitution into the equations minimizing this functional of a — qip, (3 = pip and 
9 = 0(r), $ = $(r) immediately satisfies two of them; the remaining two are as follows: 

9 sin 9 
4(9 rr + — ) = sin^ 2 H — [p 2 cos 2 $ + cos9(q + psin $) 2 ]; 



$ r 0„ ^ pCOS$. 2 . 2 . _ 

$ rr H = -ctg-0 r Q r H - 2 — (gcos --psm -sin$). 

T' Zi T' Zi Zi 



(32) 



The numbers and are to be clearly integer. The energy takes the following form: 



n2 ■ 2 0^.2 ^ 

W = 9 2 + sm 2 -$ 2 r + — 



p 2 9 

[q 2 + 2pq sin $ + p 2 ] H — - sin 4 - cos 2 $ 



(33) 



To find solutions with the finite energy one must choose the conditions 9 — > kn at the infinity 
with k being integer. The second equation implies that at the infinity $ — > (2m + l)n/2 
when 9 ^ kn. For the simplest case with minimal different values of the winding numbers: 
p = 1, g = 3 we have found the exact solution 

A 2 -r 2 



sin $ 



A 2 +r 2- 



cos 0(r) 



r 4 (Z 2 + A 2 ) - (r 2 + A 2 )/ 4 
> 4 (/ 2 + A 2 ) + (r 2 + A 2 )/ 4 ' 



(34) 
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The energy of the solution does not depend on two radii A and I and equals to 16tt (half 
the energy of the electrical skyrmion). The asymptotic behaviour of the spinor is (two cases 
correspond to two signs for cos 6) at zero and at the infinity 



/n\ / ex P( 2 ^) \ / cxp(2iy) x 



from 







to 



\ / 



, or from 



exp(i<£>) 

V o / 



to 





W 



respectively. 

Now let us discuss the solutions from other point of view. In SU(3)-isotropic case the 
density of energy is invariant under the transformations 



f = exp[z^ 7fc A«] 



k=l 



generated by {A^}, where eight parameters 7t do not depend on the coordinates. If ip2 = 
Ttpxi then 

8 8 

J2(^ (k w = E(^A (fe V 2 ) 2 

k=l k=l 

at each point. In other words, both ipi and ip2 are the solutions of variational equations with 
the same energy. Hence, all the solutions can be classified on the irreducible representations 
of SU(3) group, or supermultiplets. Generally speaking, two spinors ipi and ip2 from the same 
supermultiplet represent only one solution in different parameterizations. For instance, five 
solutions (iz) — (v) represent one solution {-01,^2,0} which is in fact a well known Belavin- 
Polyakov solution constructed on two components. 

The solutions like electric skyrmion (see Eqs.(p0|)-(p2[)) and solution (vi) from this section 
correspond to another supermultiplet and can be transformed into each other by a rotation 



around z-axis. Solutions with the R vector distribution (|23| ) derived in the beginning of this 
section can be reduced to electric skyrmion by the transformations generated by A^ 4 ) and 
A®. Another non-trivial supermultiplet can be constructed from magnetic skyrmion (^6|). 
More detailed analysis of the SU(3) supermultiplets will be given elsewhere. 

Concluding this subsection we have to notice that the dipole-quadrupole solutions can 
be obtained also in more general J\ ^ J 2 model but the corresponding equations due to 
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the term proportional to J2 — J\ cannot be solved analytically similarly the situation in 
S=l/2 model with exchange anisotropy for the out-of-plane vorticesll. It should be noted 
that varying the "anisotropy" parameter A = J2/J1 from A = to A — > 00 we deal with a 
transformation of purely magnetic solution to a purely electric one. One might expect that 
magnetic skyrmion (|T^) would be stable given < A < A Cl , while electric one would be 
stable given A C2 < A < 00. In the intermediate range A C1 < A < 1 and 1 < A < A C2 we expect 
the stability of dipole-quadrupole skyrmions with predominant "magnetic" and "electric" 
behavior at infinity, respectively. Given A = 1 we come to the fully SU(3) isotropic solution. 
The critical values A Cl 2 of the "anisotropy" parameter should be calculated numerically 
similarly to S=l/2 anisotropic Heisenberg modefi 



IV. CONCLUSION 

In terms of spin coherent states we have investigated topological defects in 2D 5=1 
(pseudo)spin quantum system with the bilinear and biquadratic isotropic exchange in the 
continuum limit. The proper Hamiltonian of the model can be written as bilinear on the 
generators of SU(3) group (Gell-Mann matrices). Knowledge of such group structure enables 
us to obtain some new exact analytical results. The analysis of the proper classical model and 
its topology allows us to get different skyrmionic solutions with finite energy and the spatial 
distribution of spin-dipole and/or spin-quadrupole moments termed as dipole, quadrupole, 
and dipole-quadrupole skyrmions, respectively. Among the latter we would like to note the 
in-plane vortices with the in-plane distribution of spin moment, varying spin length, and the 
non-trivial distribution of spin-quadrupole moments. 

One should note that for traditional spin systems like 3d magnetic oxides with S > 1 
the biquadratic exchange as a rule two orders of magnitude smaller as compared with usual 
Heisenberg bilinear isotropic exchange. It seems, for such systems the above analysis may 
play only purely theoretical interest. However, for systems with the non-quenched orbital 
moments, (pseudo)-Jahn- Teller effect, and other forms of (pseudo) degeneracy the effective 
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pseudo-spin Hamiltonian given S > 1 may include different large non-Heisenberg terms. 
Namely for such systems we may expect the manifestation of unusual topological defects, 
including those addressed above. 
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